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Chapter 1 

Introduction 



The aim of the present work is to show how the algebraic approach to the 
question of integrabihty can be given geometric foundations. The notion of 
first integrals for dynamical systems (Hamiltonian in particular) is almost 
always formulated with the tacit assumptions that the underlying phase 
space is Euclidean (or given the additional symplectic structure). The con- 
cept of the metric structure, crucial for Manifolds, is usually omitted - and 
no wonder, since there is no clear way to introduce a distinguished norm for 
a general system. 

For Lagrangian mechanics there is an equivalent description by means 
of the Jacobi or Eisenhart metric with which the flow of the system can be 
made geodesic with respect to precisely determined notion of length [ 211 [5)1 . 
However, this only takes into account the configuration space of the system, 
not the whole phase space, which means that for a Hamiltonian system one 
only considers the coordinate subspace with second order equations on it 
and nothing is said about momenta. This is unacceptable when the charac- 
teristic exponents, or chaos is to be investigated - the possibly exponential 
divergence of the trajectories has to be analysed in the full phase space. 

Instead of trying to find one canonical structure, this work presents the 
general view that differential geometry has to offer in that topic. A dy- 
namical system is analysed on a Riemannian manifold, to show how the 
well-known equations and definitions have to be modified, with a special 
view to the questions of integrabihty and Lyapunov exponents. The latter 
are examined in detail to obtain a differential equation whose solutions are 
the so called "time-dependent" exponents, which tend to the standard ones 
in infinite time. 

The integrabihty chosen for study here is generally understood to be the 
existence of enough first integrals (in involution for the Hamiltonian sys- 
tems), which are meromorphic functions over the complexified phase space. 
On the one hand, this bears clear consequences on the geometric picture 
and is easily translated into the language of differential geometry. On the 
other, the question of proving such existence leads to the deeply algebraic 
properties of the system such as the analytic continuation of the solutions 
in complex time, the differential Galois group and solvability of linear dif- 
ferential equations by quadratures. 

The group-theoretic approach relies on particular solutions and com- 
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putations of explicit equations in some coordinate systems. As mentioned 

before, it is usually silently assumed from the beginning that the mani- 
fold is Euclidean which leads to significant simplifications when it comes to 
derivations, vectors and matrices. Differential geometry requires that we 
use covariant derivatives, tensors and distinguish between 1-forms and vec- 
tors. It seems like an unnecessary complication to add but as it turns out, 
many theorems are much more straightforwEird to prove (like the existence 
of integrals of higher variational equations for example), some concepts like 
self-adjoint operators appear naturally and allow for applications of known 
theorems from other branches of mathematics £md last but not least, it is 
possible to identify which structures are in fact of geometric origin, which 
can be defined as coordinate invariant and which are only justified by the 
efficiency of calculation. 

The main body of this thesis is divided into two parts. In Chapter 2 the 
differential Galois group fundamentals are explained, and the basic defini- 
tions and steps in investigating integrability are described. Next, Chapter 
3 deals with the geometric approach, showing how the concepts introduced 
earlier need to be changed, how general dynamical systems are described 
in this language and finally showing how Lyapunov exponents can be con- 
sistently introduced in a covariant way and how this leads directly into an 
easily applicable numerical routine. Since the notation is index-free (stress- 
ing the independence from coordinate systems or even coordinate bases) al- 
most all derivations are given in detail. The reason for this is first to make 
it possible to follow exactly the flow of exposition if the reader so wishes, 
and second the fact that there seem to be very few practical applications of 
the notation. Hopefully, the successful formulation of the basics of phase 
space dynamics presented here will prove it is not reserved for pure math- 
ematics only. 

The third part of this work contains three examples on applying the 
Lyapunov exponents formula or algorithm, and simple comparison of the 
algebraic formulation of the normal variational equations versus the same 
equations as obtained for the simplest Euclidean manifold with the Levi- 
Civita connection. 
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Chapter 2 

Algebraic Approach to 
Integrability 



Integrability is mostly synonymous with existence of first integrals - func- 
tions that are constant along the solution. However, for some classes of 
dynamical systems there are additional, specific requirements. 
Consider a general autonomous system 

Xi = ^=Viix), i = l,...,N. (2.1) 

as 

It is said to be integrable in the Euler sense if there exist A^— 1 first integrals 
Ji, which mesms that for i = 1, . . . , N — 1 

3 T 

Ji = -^Vk = 0. (2.2) 

OXk 

Note that only N — 1 integrals are required, since they formally reduce 
the system to a one dimensional one, which can be solved by a quadrature. 
That last step gives rise to another constant, which can be considered as 
the A'^-th first integral. In the case of autonomous systems this reflects the 
freedom of translating the solution "in time", or along the trajectory. 

Most physical systems posses additional structure of being Hamiltonian. 
The dimension is then necessarily even N = 2K, with the first K coordi- 
nates customarily denoted q, , and the other K — called momenta — denoted 
Pi. The system then has the following form 



dH 

Qi = 

dpi 

Pi = -^, i = i-,...,K, 

OQi 



(2.3) 



where H is called the Hamiltonian of the system. 

While the first integrals are defined exactly as above, these systems are 
such a restricted class, that much less is needed for integrability. We say 
that s Hamiltonism system is integrable in the Liouville sense when it has K 
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first integrals which are in involution. The additional requirement means 
that for any two Ji and Jk 



[J- Jk] — y^( — — - ] - (2 4) 

" ^ \ dqj dpj dpj dqj I 



By construction the Hamiltonian is itself a first integral, and is in involu- 
tion with any other additional first integral, as 

[J, H] = j. (2.5) 

These are only the basic facts and notation needed here, and a com- 
plete exposition of the topic can be found in |[T1[6]|. The aspect of algebraic 
formulation of integrability that is of main interest in this work is the ex- 
istence of first integrals and the criteria or tests of this property. There are 
no general algorithms for finding constants of motion explicitly, which goes 
hand in hand with the fact that solvable systems are a rare exception (zero 
measure set) among all dynamical systems. There are however conditions 
which integrable systems must satisfy (i.e. necessary conditions) and that 
allows of excluding most systems so that only a few particular cases poten- 
tially solvable are left. They can then be subject to a more detailed analysis 
not possible for a whole general class of systems. The method described in 
the next section deals with necessary conditions of existence of first inte- 
grals and is one of the most restrictive. It has been successfully applied to 
many systems, determining completely the cases which are non-integrable. 
For example, it has been proven that a Hamiltonian system with a homo- 
geneous polynomial potential and natural kinetic part of any dimension is 
only integrable in at most finite number of cases II17II . For other examples 
of application see lfl6l[l9l[T5J . and Ii14l i2iI for a detailed introduction. 



2.1 Differential Galois theory 

The first step required for this approach is to linearise the general equa- 
tions ( I2.lt around a particular, non-constant solution ^(s) so that 

x,{s) = + exf^ + ^ + . . . , (2.6) 

where e is some small parameter, and each is called the n-th variation. 
The original equation then yields in each order of e the n-th variational 
equation 

■^"'^ (2.7) 

with /(^) identically zero, and in general polynomial in the variations. It 
is the first variational equation (VE) that will be of most interest, although 
the higher ones are also influenced by integrability providing further nec- 
essary conditions |'14'|. 

As the VE is linear it has N linearly independent solutions, and one of 
them is simply t/j. If one imagines that the components of x'^' represent 
separation of two nearby trajectories of the full system then this trivial 



■ (ri) ^ 
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solution corresponds to the displacement relative to the same trajectory 
translated in the independent variable s. It is called the tangential part of 
the VE. Although it is easily solvable on its own, the complete VE might 
still as a whole indicate non-integrability. However, there are N of them, 
and for the sake of practicality, one is almost often forced to work with the 
components left after discarding the tangential part. This fact gives rise to 
the normal variational equations. 

Formally, the normal bundle is defined to be N^,7\/ ;= T^JM/Tijj, where 
for brevity ip also denotes the trajectory associated with the particular solu- 
tion tp{s), and M is the base manifold of the original system. The projection 
TT : T^M ^ N^M is then used to construct the normal variational equations 

y = ^ (^1^ ° TT-i (y)^ , y e N^M. (2.8) 

In this approach, there is no metric structure so, strictly speaking, there is 
no orthogonality of the normal components to the tangent ones. Usually, 
the particular solution appears when the first iV — 1 dependent variables x 
are equal to some constants and ip is determined by 1 differential equation. 
Assuming the Euclidean structure (of both the base manifold AI and thus 
also of the fibre of TAI), the tangential part of VE is associated with the 
direction x n and the remaining directions are treated as normal. 

With Hamiltonian systems, the reduction involves one more step, as 
there is the first integral H, which can also be used. One more degree 
of freedom is eliminated by considering the system i2.8i on the constant 
energy hyper-surface. In terms of the linear variations this means 

diJo7r"^(?;) =0. (2.9) 

This time the particular solution is sought for so that N — 2 variables are 
(e.g.) zero and a pair of andpK are left to provide a second order equation 
for ^p. The above reduction by 2 degrees of freedom then gives the NVE 
involving only qi and pi for i = 1, . . . ,K — 1 (where, as before, N — 2K). 

In any case what is left is a set of (non-autonomous) linear differential 
equations 

y = A{s)y, (2.10) 

for which it is possible to define the monodromy group M . When the inde- 
pendent variable is considered as complex, one can ask how the fundamen- 
tal matrix of solution changes after analytic continuation of the solutions 
in closed loops around a point sq- Since it must be a linear function of the 
initial fundamental matrix, a matrix multiplier is obtained, and after con- 
sidering all loops, one ends up with a whole group. This subgroup of GL(C) 
will be an image of the fundamental group tti{^, so) of the Riemann surface 
^ defined by the solution ip. 

The existence of a first integrals J of the main system implies there is 
also a first integrals of the (normal) variational equations 

^ = const, (2.11) 

provided the gradient does not vanish on the trajectory ^p, but even then 
the higher derivatives yield a first integral (they cannot all vanish for then 
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J would trivially be zero). This will be presented in more detail in the next 
chapter, as at present it is enough to notice that the above formula gives a 
function / such that 

= (2.12) 

for all ginM. Such / is called a first integral of the group and it was shown 
by Ziglin | 23 1 that if the main system has a meromorphic first integral then 
the monodromy group has a rational first integral. That fact alone can be 
used for the study of integrability, but since the above can be refined still, 
let us pass to the Galois theory itself 

For a linear system 02.1OD . with the coefficients in some differential field 
K 3 A{s) (for example the rational functions of s, C(s), with the standard 
derivation of d/ds), the solution almost never lies in K but in a larger field 
F D K. It the field extension is generated by all the linearly independent 
solutions of the given equation, the field F is called a Piccard-Vessiot ex- 
tension of K. It is (up to an isomorphism) unique, and allows to define the 
differential Galois group Q of O2.10l >. as the group of automorphisms of F 
that leave elements of K fixed and commute with the derivation. 

The Galois group is bigger than the monodromy group, and in the special 
case of the equations being Fuchsian, M is dense inQ. It is still an algebraic 
subgroup of GIj{K'), where K' is the field of constants of K (C in most 
applications). It is still the case that there exist integrals of the Galois 
group, when there are integrals of the dynamical system. This leads to two 
important properties. 

1. When the system is integrable in the Euler sense with meromorphic 
first integrals, the Galois group of the normal variational equations is 
finite. 

2. For Hamiltonian systems which are integrable in the Liouville sense 
with meromorphic first integrals in involution, the identity compo- 
nent of the Galois group of the (normal) variational equations is abelian. 

The above are fundamental theorems of the theory, but for practical ap- 
plications one is usually interested in the consequences for the solutions of 
the NVE. In the former case the solutions lie in an algebraic extension of 
the field K, and in the latter in a (generalised) Liouvillian extension. Such 
an extension is a formally defined concept of a solution in a "closed form" 
and arises from K after finitely many steps each of which consists of includ- 
ing one of three types of a new element: one that is algebraic over K, its 
derivative lies in K or the derivative of its logarithm lies in K. 

In physical applications, the NVE are usually of the order two, and there 
is a general tool for checking the above property: the Kovacic algorithm 
|[Tp1. There are also particular results for special families of equations like 
the Riemann P-equation |8| or Lame equation |4|, which give explicit condi- 
tions on the equations' parameters for their solution to be Liouvillian. The 
examples cited in the previous section show how effective the method is - 
reducing a general problem of the existence of first integrals to the question 
of solvability of linear differential equations. 

Before the other (geometric) description is presented, one remark is in 
order. Although the Galois group approach is applied to physical systems, 
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for which one naturally considers the variables to be real, the mathematical 
"workshop" is located in the complex domain. This means that not only 
the dependent variables but also the time itself can be imaginary, and this 
has several implications. First of all, for all physical systems where one 
assumes sufficient smoothness, the lack of complex integrals implies the 
lack of real ones (as a special case). However, there are systems for which 
there exist real smooth integrals, that are not even real-analytic |7|. Thus, 
integrability is always depends on the domain considered, and the family 
of functions considered "good enough" to be first integrals. 

Secondly, some systems although non-integrable behave chaotically in 
the imaginary phase space and not in the real one. A particular example 
of that is the Gross-Neveu system analysed in [12|, which will also be stud- 
ied here with the help of Lyapunov exponents. Since the algebraic theory 
is usually studied without much input of differential geometry (as there is 
usually no need for metric structure) the exact connection between inte- 
grability and geometric concepts like exponential separation of trajectories 
or topological chaos remains elusive. The next chapter is thus devoted to 
treating dynamical system with more detailed differential geometry and, as 
an example, giving a new formulation of the Lyapunov exponents together 
with some implications that first integrals enforce and practical results. 
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Chapter 3 

Geometric View 



Throughout this chapter the index-free, operator notation will be used, and 
the convention is mostly that of [ 18 1 or 1 9 1 - the reader can find in depth 
introduction to the subject in those two books. 

Vectors and 1-forms will be denoted by bold symbols, and A, B, X, Xi, 
X2 and so on will mean arbitrary vectors and cf) will signify any 1-form 
(usually used to define an operator or prove an equality). Vectors will be 
considered to act on functions as differential operators so that in a coordi- 
nate map u 

i—l 

where N is the dimension of the base manifold. 

The covariant derivative V is defined for (^) type tensors (the exterior 
tangent bundle strictly speaking) by the property 

W(j3(E)X) = df3(E)X + {^1)''(3AWX, (3.2) 

with (3 being a fc-form (or a function for k = 0). And the connection itself is 
fully characterised by its action on any basis {Ei} 

yEi^J2^v®Ej^ (3.3) 

3 

where is a matrix of 1-forms (not a tensor). It allows to define the torsion 
tensor, the curvature endomorphism and the Ricci tensor as 



T(Xi,X2) 
R{Xi,X2)A 
n{A,B) 



(VxiVx2 - Vx.Vxi - V[xi.x.])^ (3.4) 
tr(X ^ R{A,X)B). 

With the introduction of the metric structure (a symmetric tensor) g 

0(Xi,X2)^ (Xi,X2) eM, (3.5) 

the Levi-Civita (torsionless and Riemannian) connection can be introduced 
by taking 

T = 0, 

d (Xi, X2> = (VXi, X2> + (Xi, VX2) . ^^'^^ 
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The second requirement implies also 

d77(Xi, . . . , Xjv) = ^ viXi, . . . , VXj, . . . , Xn), (3.7) 

j 

for the volume form 77. The metric also allows to define the Riemann tensor 

R{Xi,X2,A,B) - {R{Xi,X2)A,B) . (3.8) 

The interior product t is defined as 

{LxO(Xi,...,Xk) ■.= C{X,Xi,...,Xk), (3.9) 

where C is some (fc + l)-form. The following property of the Lie derivative 
£ will also be useful at some point 

»exC = d(txC) + ^x(dC). (3.10) 

3.1 The variational equations 

Consider now the following construction (Figure l3?T] ). Given the field V and 
a curve 7(0, 1) (not being one of the field's integral curves), we can construct 
the image of ^{l) under the flow of V with s = si. The image curve will be 
denoted by j{si,l), and by using the same / we understand, that each point 
of the new curve is the image of a point of the original curve for that value 
of /. Thus, the parameter /, which needs not be the natural parameter, gives 
rise to a new vector field Z ^ tangent to family of transformed curves 
(along with the original one). 




Figure 3.1: A visualisation of the family of curves 7(5, 1) generated by the 
fields V and Z (with their dependence on the point indicated). 
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By construction, for any two points p and q on 7(0, 1) carried into r and x 

on 7(si, the differences of their I values are the same. Or, in other words, 
that the flows generated by V and Z (g and h respectively) commute 

g,, {hi, ip)) = hi, {g,, ip)) ^ [V, Z] = 0. (3.11) 

Taking into account that the connection is torsion-free, this will allow us to 
transform the derivative of Z along 7(5). 

T(V, Z) - [V, Z] = VvZ - VzV = ^ ^Z := VvZ = VzV, (3.12) 

OS 

which is, in fact, the first variational equation (VE) 

Z = (Vy)(Z), (3.13) 

where the dot will, from now on, denote the covariEuit derivative along the 
field V. 



3.2 Projected variational and deviation equa- 
tions 

Let us introduce the operator h projecting a vector field on the normal (with 
respect to V) bundle NM 

h = l-v^V, (3.14) 



where i; is a one-form dual to V given by 

{V,V)' 



v{X) = i^^, (3.15) 



so that v{V) = 1. Obviously we have 

/i^ = h, h{V) = 0, (3.16) 

and for constant norm of V 

{VV,V) = =^ v{S/V) =0. (3.17) 

The norm will also be written shortly as l/A/" := {V,V), and will be included 
explicitly in calculations so that all the formulae hold also in the Lorentzian 
case M = —1. For brevity, a projected vector will also be denoted by 

_lZ h{Z). 

Consider now the derivative of the projection along the curve 7(5) 

^^i_Z = Vv{Z-v{Z)V) 

= WzV-v{Z)VvV-{Wvv{Z))V (3.18) 
= V(z_^,(z)v)V - {Wvv{Z))V 
= V^zV -V{v{Z))V, 
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where we have used the commutation property ( I3.12t . Acting with h on 
both sides of the above equahty we obtain 

±^±Z = A^^zV. (3.19) 

OS 

To obtain the deviation equation we differentiate again 

= R{V, i_Z)V + V^zV + V[v,,z] V, 

and project onto NAf 

±^±^±Z = h {R{V, i_Z)V) + A^^zV + ±V[v,^z] V 
OS OS I J 



(3.20) 



= RiV, i_Z)V - N{V, R{V, xZ)V) + ^ ^zV + ^yv.^z\V 
= Riy, i_Z)V + aV^zV + ±V[v.^z] V. 

(3.21) 

The commutator can be further simphfied 

[V,^Z] = [V,Z-v{Z)V] 
= [V,-v{Z)V] 
= -V {v{Z))V -v{Z)[V,V] 
= -J\fV{(V,Z))V (3.22) 

= -M ( (v, z) + (V, VzV)) V 



= -Af{V,i_ZjV, 
the scalar coefficient above is 

V{v{Z)) = U{(^V, + (v, Z^) - {V, Z) 

= N (iy, i_z) + (V, v) v{Z) + {V, VzF) - 2 (v, V) v{Z) 
= N ((y, i_z) - (V, v) v{Z) + {V, V^zV) + v{Z) Iv, V 

^uUv,i_z) + {v,v^zV)). 



(3.23) 

Substituting the above, the second derivative becomes 

= R{V, i_Z)V + ^y^zV - U (iV, i_z) ^V. (3.24) 

This can be written in a shorter form, which allows further analysis 
to be simpler, upon introducing the Fermi derivative. Since V does not 
necessarily define a geodesic flow, we want to define a new derivation along 
7(s) which would satisfy the following properties: 

Dp 

1. -Tfv = Q. 

OS 
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2. = ^X, for geodesic V. 

OS OS 

Dp Dp d 

3. If— ^X = = -^(7,then^(X,C/) =0. 

OS OS ds 

4. ^^X = ±^X, for X orthogonal to V. 

OS OS 

The last property allows us to find the explicit form of the Fermi derivative 

^X ^ ^X -U(V,X)V +M (v,x)v. (3.25) 

OS OS \ / 

Equations 03.19l > and ( I3.24D now read 

Dp 

-jfi_Z = ^V^zV =: A_AUZ), 

A (3.26) 

-^^_Z = R{V., i_Z)V + ^V^zV - N l^i_Z., ^v) =: 

These can be called the projected variational and deviation equation respec- 
tively. For both operators introduced here we have 

aA, $ : TM ^ NM, (3.27) 

where the normal bundle is taken with respect to V , so that the matrix 
of components of \A explicitly only has non-zero elements for directions 
orthogonal to V, and is effectively a (A^— 1) x (A^— 1) matrix (when acting 
on NA/). 



3.3 The Raychaudhuri equation 

Usually one decomposes A into its symmetric and antisymmetric parts us- 
ing the corresponding index notation 

^m;- = 1^(^;.) + V[^;H- (3.28) 

However, with the assumed definition of the covariant derivative, A is not 
a tensor of type (2), but can rather be identified with a tensor of type (J). 
The question of "transposing" such an operator can be dealt with naturally, 
when one recalls that the metric tensor defines the musical isomorphism, 
so that the following diagram commutes 

A — ^ \>A 

1 1 It , (3.29) 

A^ {\)AY 

where the new symbols above are defined as 

{\>A){X,U) = {A{X),U), 

{ct>'^)(X,U)^4>iU,X), (3.30) 
((tt</.)(X),(7> =0(X,(7), 
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It is, of course, possible to only define the Hermitian adjoint operator of A 
as 

{A^X),U) := {X,A{U)), (3.31) 

but then the analogy with simple transposition is not as clearly visible. 
As can be seen the adjoint is really a transposition of arguments in the 2- 
form (j) canonically isomorphic to the given endomorphism A. In the index 
notation the above operations take the form 



(3.32) 



The adjoint is especially simple for operators of the form X as can 
be directly checked 



= ci>{B) {A,X) 

= \>X{A) {U,B) 

= {{\>X®U){A),B). 



(3.33) 



or(08)X)t = (bX)®(fl0). 

We are now ready to decompose _lW into its self-adjoint (Hermitian) 
and anti-self-adjoint (anti-Hermitian) parts. This is essentially the sym- 
metric splitting of the associated 2-form (j). The reason for carrying out 
this procedure in the operator approach is that we are dealing with lin- 
ear differential equations and any eigenvalue problem will be much more 
straightforward. We will be able to use all the standard theorems regarding 
Hermitian operators. 

Decomposing now with regard to f we have 

^^Z = e{^Z)+u:{^Z), (3.34) 

with 

= iA" ■.= \{aA + aA^), 

f (3.35) 

so that 

6it = 61, J = -CO. (3.36) 

One could think of 9 and to as generating the flow's expansion and rotation 

respectively. 

We turn next to the deviation equation to obtain the derivative of ±A 
§.Z = *(.Z) = 5f (^,Z,) =: 5|i4,.Z) + ^ (|£.Z) , <3.37) 



which we rewrite as 

Df±A Dpd Dpuj 



<^-±A\ 



ds ds ds ' (3.38) 

= $-6»2-a;2-{6(,a;}, 
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where {,} stand for the anti-commutator. It is straightforward to con- 
struct the Hermitian decomposition of the left-hand side of the above equa- 
tion, but one must ask if the derivatives of 9 and uj are themselves the 
(anti-)Hermitian parts of the left-hand side. Or, in other words, if 



ds \ ds J 



(3.39) 



First we prove this holds for the covariant derivative (which is expected for 
a Riemannian connection) 



t 

\ ds J ' / \ ^ ds 



V{{X,A{U))) - (^,^(C/)) - {^'(X),^ 
Vi{AHX),U))-(^A^ (^)'^ 
-V{{A\X),U)) + (^AHx),U 



s 



(3.40) 

For the Fermi derivative it thus suffices to check the behaviour of the last 
two terms, which are all of the form 03.331 1. They can be written as 

Wi = -Af(\>V (E)V -bV ^V), 

so one immediately has w} = —W\. The question ( I3.39D can now be rephrased 
as 

iy^^o A- AoWx)^ = Wio A^ - A^ oWx. (3.41) 
As {A o _B)^ — b'^ o , the above is identically satisfied for this particular 

As we assume zero torsion, and the connection to be Riemannian, we 
have for the Riemann tensor 

R{A, B, X, Y) = R{X, Y, A, B) 

so that using an auxiliary operator 

Wr{X) = R{V,X)V 

it is possible to obtain the adjoint of i? as a function of its second argument 

(WRiX), U) = R{V, U, VX) := {R{V, U)V, X) 

= {X,R{V,U)V) (3.42) 
^{X,wI{U)) = {X,Wr{U)), 
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and it turns out to be self-adjoint. This allows us to write 

= Wr + {aVV^ - MaV ® aV, (3.43) 

or 



$^ Wb. + (±VF)^ - J\f{\>±V) (8) ±V, 

$^ = (±vy)^. 

Together with equation 03.381 1 this gives 



(3.44) 



2 , ,2 



ds 

iAVV)^~{9,Lo} 



(3.45) 



ds 

We introduce here a new quantity -d ~\,xd 

tr^ = ^(tr0)=,?, (3.46) 
as as 

since the Fermi derivative of a scalar function is just the derivative with 
respect to s. Taking the trace of the right-hand side of the first of equations 
(IHSS) yields 

^T^i-xr T/•^ , ^-„/ V7T>^ Kf Ur tV\ ^^.lai , . .2^ 



= n{v, V) + tviyv - M (vv, v'^®v) -m(v,v'^ - tT{9^ + u^) 
= n{v, V) + tr(VF) + TV (y, Vvv'^ - N (v, - tv{e^ + u?) 

= 7^(y, V) + tr(Vy) - tr(6'2 + u^)^ 

(3.47) 

where we have used the fact that 



V {v,^! = 

and that 

tr(</.(g)X) ^ (j}{X). 

Writing conventionally 

^ -1 



N-l 

to separate the trace-free part a, one finally arrives at the Raychaudhuri 
equation 

^ TZ(V, V) - tr - tr ct^ !_ij2 _^ tr(VF). (3.48) 

ds A* — 1 

3.4 Higher order variational equations 

Let us turn now to the flows generated by the fields V and Z as shown 
on Figure 13.11 Let u : AI —> denote a local set of coordinates, and 
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consider a function / at point r = gs{p), where the locEil diffeomorphism 
is associated to the vector field V so that 



=: Vr{u) = fr. 



(3.49) 



These are simply the field's components in the coordinate basis associated 
with u (not that the subscript indicates the point, not a vector component). 
Taking their derivative with respect to s we write 

"^-fr = Vrif), 



or 



— u{g,{p)) = Vr {Vr{u)). 



SO that 



d" 
ds" 



(3.50) 
(3.51) 
(3.52) 



Take now another point q — hi(j>), generated with the flow of Z, and expand 
its coordinates, as is always possible in a local map 



u{q) = u{hi{p)) = 



n=0 



n! 



dZ' 



1=0 



n=0 

oo 



(3.53) 



n=0 



We define the sequence of displacements w„ to be Z^^u), and call it varia- 
tions of the n-th order. The partial sums correspond to points p„ such that 
uiPn) = Z]i<„ Ui, and the n-th variation is then a vector in (coordinate 
space) connecting the point p„ i to 

Thus, only one vector field is needed to describe all the variations, al- 
though we can formally write the higher order equation (HVE) as 



^t.„ = y(z(")(u)) = z(") {v{u)), 

by Emalogy to the originsd dynamical system 

u = V{u) <^ u^' = V^'. 
For example, the first VE equation is 



(3.54) 



ds 



(ZK)) = Z{V{un) = Z{Vn, 



Z'' = Z" 



(3.55) 



and the second 



-{ziz{un)) = z{z{vn), 



= z^z' 



Z" 



dZ^ dV 



(3.56) 
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or, in the usual notation, 



\ o (3-57) 

The whole construction is evidently coordinate-dependent, as the chart 
u enters explicitly into the definitions, and although the first variation can 
be made into a single vector Z = u'^d^, the same cannot be done for the 
higher variations. It follows from the fact that a single field Z is enough to 
define the diffeomorphism hi which gives the full transformation of one tra- 
jectory onto another, for a small but finite separation /. In practise however, 
we only know the values of Z on the particular trajectory 7, and not its de- 
pendence on the coordinates u, and we cannot construct the derivatives of 
Z (except in the V direction). That is why the Taylor series of 03.53D is anal- 
ysed term by term, each with its own equation. Note that a HVE of a given 
order n requires the knowledge of the solutions of all the lower HVE's up to 
n—1, which makes the equations non-homogeneous (and non-autonomous). 



3.5 First integrals 

With the present notation, a first integral of the system is a function J : 
M -^R such that 

^ J = V(J) = dJ(V) = 0. (3.58) 
as 

A key fact regarding these quantities is that, if there exists a first inte- 
gral J of the original system then there also exist first integrals J„ of the 
variational equations of all orders 

J„:=Z(")(J). (3.59) 

To see that these are in fact first integrals of the VE, we make use of the 
commutation property again 

= y(z(")(j„)) = = 0. (3.60) 

as 

Just as the HVE, the integral explicitly involves the solutions of the HVE's 
of lower order. 

However, it can so happen, when taking a particular trajectory 7(5) of 
V along which the variations are considered, that 

dJ|^(s) EE 

or in general that 

aV|^(^)=0, for j<m, (3.61) 

where m must necessarily be finite or the whole first integral would be a 
constant function, and is taken to be the smallest integer with the above 
property. Note first that d"\J is a well defined tensor of type (^) because 
all the partial derivatives of orders lower than m vanish so that only the 
highest derivative is left 

(5™ J)^(,)(Z__^) = Z{... Z{J) . . .)|^(,). (3.62) 

m m 
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The reason to choose d to denote this derivative is that V is understood to 
act on the (possible exterior) tangent bundle, while for d one has = 0. It 
is to be remembered, though, that we do not associate d with derivation in 
any particular coordinate system, as can be seen in ( I3.62l >. 

Let us take now a VE of n-th order with 5" J = 0, and a following func- 
tion of the n-th variation 

Z(")Z("''J = Z(")j„/ (3.63) 

where n' is such that n-\-n' = m constructed as above. Then {d'^Jn')'~,{s) is a 
first integral, expressed explicitly as a function of the n-th variation. 

Imagine now, that the system has N functionally independent first inte- 
grals in a neighbourhood of 7(5). This means that there are N independent 
first integrals of the first variational equation (insert Ziglin Lemma), and 
if the original integrals are independent on 7(5) itself, this means they all 
have non- vanishing associated 1-forms dJ^. Consequently we have N inde- 
pendent vectors 

= tt(dJO, (3.64) 
with which the integrals of the variational equation can be written as 

J,4 = Z{ J,) = dJ,;(Z) = (r„ Z) . (3.65) 

Since that expression is constant 

Q = {VvY,,Z) + {Y,,VvZ) 

= (r„z) + (r„VzV) (3.66) 

= (y„z) + ((vy)tF„z), 

or 

^Y, = -{VV)^Y,, (3.67) 
OS 

and we say that Y^ satisfied the adjoint equation to the variational equa- 
tion. In fact, we have just shown that N first integrals give us a basis of 
solutions of the adjoint equation. From them, the solution of the VE can be 
obtained by means of the fundamental matrix (or rather operator). 

Let F be the fundamental operator of the VE and the adjoint equation 
respectively, defined as follows 

N 

F = ^Z*o(^Z,, (3.68) 

i=l 

where ZiQ are the initial conditions of a basis of solutions {Zi}, and {Z*} 
is the dual basis. The above operator acts on a constant vector to yield 
a solution of the variational equation with the constant vector as the ini- 
tial value. When the constant vector is prolonged (along the trajectory 7) 
with the equation Vv^io ~ (which is not the same as the VE), then the 
derivative of the above operator is simply 

^F = AF. (3.69) 

OS 
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Another way of making the initial values global is to require that the field 
{Zi} is, at some point, equal to a basis {-E,} which is globally parallel: 
VEi = 0. This is the usual case with the implicit assumption of the base 
manifold (and connection) being Euclidean. 

Define now P to be the fundamental operator of the adjoint equation, 
and Yio be the initial conditions of the appropriate vector fields; then 

c,, = {Y,,Z,) = {FY,o,PZ,o) 

^{P^FY,o,Z,o)^{Y,o,Z,o), 

where the last equality is the consequence of this scalar product being con- 
served as the first integral. As the initial conditions are arbitrary, this 
means that 

P^{F^y\ (3.71) 



3.6 Normal variational equations 

As was mentioned in the chapter on algebraic theory, the variational equa- 
tion can be reduced in order by 2, when a first integral is known. We have 
already seen how the first part of this reduction works - by projecting the 
variational equation on the subspace orthogonal to the trajectory (tangent 
vector). The second step is carried out almost identically, only this time the 
vector used for projecting is 

Y = ttdJ, (3.72) 

where J is the known constant of motion. By definition Y is orthogonal to 
V 

{Y,V) ^dJ{V) ^0. (3.73) 
Thus ±_Z can be further decomposed as 

±z = ^Jy^y) ^ + ^^-^^^ 

The last vector _u_Z is tangent to the hyper-surface of constant J and orthog- 
onal to V . Just as Z{J) is the first integral of variational equations iZ{J) 
is a first integral of the projected equations because ( I3.73D implies 

{Y,Z) = {Y,^Z). (3.75) 

The reduced variational equation can be obtained from ( I3.74D by tak- 
ing the covariant derivative Vv and projecting on the subspace orthogonal 
to both V and Y or, which amounts to the same, by applying the Fermi 
derivative 



= ^AUZ) (2 {a_aHy),y) Y {Y, Y) ^AHY)) , 



(3.76) 



and projecting with respect to Y 



l^Z . ^A L Z + ^^^y) + ^AHY)^^, (3.77) 
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where the Fermi derivative of Y is known from 03.67D . so that finally one 
gets a familiar looking equation 



il^alZ = alAUZ) + 2^9{Y)^^^^. (3.78) 
OS [Y ,Yj 

In contrast with the (once) projected equation, the above contains an addi- 
tional term, which means it is coupled with the degree of freedom parallel 
to Y, albeit the "coupling" (Y, j_Z) is constant, as it is the first integral. 

When it comes to the normal variational equation, one is interested in 
the full derivative of the variations, not only its transverse value, so that it 
is necessary to go back to equation ( 13.181 1 and change it to 

^i_Z = AUZ)-AfV{{V,Z))V 
OS 



A{i_Z)-Mi^(V,Zj + (V,VzV))V (3.79) 
AUZ)-Af(v,i_z) V. 



The normal part is by definition the variation tangent to the hyper-surface 
of the first integral, so that (Y, ±Z) — 0, and one can simply write the 
normal variational equation as 

^^Ai_Z ^AUZ)~Af{v,ALz)v. (3.80) 

3.7 Lyapunov exponents 

Define new operators L and A 

L = i^tf = e2(^+i)^(^), (3.81) 

so that 



A = ^^^°g^. (3.82) 

2(s + l) s + l 

The Lyapunov exponents are then defined as the eigenvalues of 

Ao = lim A. (3.83) 

s — ^oo 

By construction L is self-adjoint, and has positive eigenvalues, so that A is 
well defined. Clearly, the operator Ao depends on the particular solution 
one uses for s tending to infinity. Or, in terms of F, it depends on the points 
in the phase space of the variational flow. There is also a slight difference 
from the usual notation, where s is present instead of s + 1. This change is 
introduced to make A well defined for s = by formula 03.811 1. 

The above formulae require the knowledge of the fundamental matrix 
F in order to be able to determine the values of the exponents. In practise 
one uses various algorithms Il22l |3l to reconstruct the spectrum, and they 
mostly rely on integration of the variational equation and, if the system is 
nonlinear, the original equation as well. 
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The following considerations provide a differential equation for the Lya- 
punov spectrum, or, strictly speaking, an equation determining a operator 
whose eigenvalues are the same as those of A and accordingly tend to the 
Lyapunov exponents with s ^ oo. This is achieved through a similarity 
transformation 

FLF-' = FF^ = e2(.+i)f Af-^ g2(s+i)£^ (3 34) 

Let us turn next to the formula for the derivative of matrix (or operator) 
exponential 

g-Afd^M^ /' e-°^Me"^^da. (3.85) 
ds Jo 

It is directly applicable to the covariant derivative (or Fermi derivative) 
thanks to formula ( I3.37l >. which could be rewritten as 



ds 



Dp 

OS 



(3.86) 



with each side of the above understood as acting on some vector. The stan- 
dard derivation now gives 



da \ ds 



ds 



which, when integrated over a leads to 03. 851 *. 
Applying the above for M = 2(s + 1)£ 



{FF^y\FF^ + FF^) = C e-2"("+i)^(2/: + 2{s + l)£)e2"(^+i)^da, 

Jo 

F-^F-\AFF^ + Fi^Ut) = 2C + 2 f e-2"('*+i)^(s + l)£e2"('*+i)^da, 

Jo 



-i)Ad£^ + y^t = 2/: + 2(.s + 1) r e-^'^^'+^^^'^^Cda, 

Jo 



(3.88) 

where the adjoint (not to be confused with the other adjoint, which will be 
explicitly denoted as "Hermitian adjoint" throughout this section) operator 
has been introduced as 

{AdW){U) U] (3.89) 

and satisfies 

e^Ue-^' = e^'^^U. (3.90) 

The integral can now be evaluated symbolically (treating the exponent as 
a scalar), which is justified by the direct computation of the exponential 
series. It gives 

„-2(s+l)Ad£ _ -ti 

(e'2(.+i)Ad£_j)^ + A + A^ = 2C- — -C, 

2Ad£ _ _J_ ^^-^^^ 
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or, upon defining 
we finally obtain 



(3.92) 



t = V(2(s + l)Ad£)(6') + \A,0\ £. (3.93) 

s + 1 s + 1 

Equivalently, by including the (s + 1) factor into the operator £ = (s + 1)£, 
we have 

£ = ?A(2Ad£)(e') + [A'C]. (3.94) 

This equation appears to be much simpler than the preceding one, but one 
has to remember that £ diverges at infinity for non-zero Lyapunov expo- 
nents. 

The above derivation requires a few crucial remarks. First of all, the 
adjoint operator, trivially has a zero eigenvalue, as any operator commutes 
with itself Ad£ is thus not invertible. Thus, it is important to remem- 
ber that the fraction notation is meant as symbolic - the "division" by an 
operator is used to write the infinite series obtained by integration as a 
simple function. On the other hand, the behaviour of the fraction at "zero" 
(a non-invertible operator) is regular because the function Xjil^ix) has a re- 
movable singularity at a; = and is taken to be equal to its limit there: 
l/?/;(0) — ?/!(0) — 1. This also makes clear, that neither V' nor l/V-! is zero on 
the real axis. 

Although l/V", when viewed as a series, has infinite radius of conver- 
gence, the same is not true for the other fraction, which is represented by 
t/i (the radius at zero is 27r). The function is, however, well defined for all 
real values of the argument because its singularities lie on the imaginary 
axis. That is why one first replaces the well-behaved series by 1 /i/i and then 
consequently uses only ^. 

Finally, the adjoint operator is also Hermitian as £ is, by definition, 
Hermitian. To see this, let us define a natural metric in the vector space of 
operators 

{A,B) =\.x{A^B), (3.95) 

so that 

tr(AUd£S) = tr(A^£S - A^BC) 

= tr((£A)tB-(^£)tS) 

' ' ' (3.96) 

= tr([£,A]ti?) 
= {AdCA, B) , 

which means that Ad£ has real eigenvalues and V is well behaved on its 
spectrum. 

Equation 03.93D is thus a matrix differential equation with the initial 
condition £(0) = 1, and involving an operator on the space of operators 
i/;(Ad£). Since Ad is linear, the adjoint can be considered as an N"^ x N"^ 
matrix acting on an N"^ x 1 vector (representing aN x N matrix). The equa- 
tion is solved for the derivative which means it is easily implementable nu- 
merically. It only requires the knowledge of a particular solution 7 around 
which the linear approximation is considered. And although the operator 



26 



C is defined with the use of the fundamental matrix, there is no need of 
obtaining the basis of hnear solutions to solve for A because their spectra 
are the same. 

The dimension of the equation seems to complicate matters a lot be- 
cause, for example, a Hamiltonian system of two spatial degrees of free- 
dom, which requires a four dimensional phase space, gives rise to a 16 x 16 
adjoint matrix. Evaluating i/; on such a matrix cannot be achieved by a se- 
ries, as mentioned earlier, and requires an eigenvalue decomposition, which 
would make the calculations cumbersome. Fortunately, this is not a general 
TV-^ X iV^ operator, and the knowledge of £ is all we need. 

Take any Hermitian operator A with an orthonormal basis of eigenvec- 
tors {Ci} such that A{U i) ~ \iU i, another operator B, constructed from 
those eigenvectors, and forms of the dual basis 

B = Ul®Uu (3.97) 
for given k and I. The action of the adjoint of A on i? is as follows 

{BA^ =AB^ ^U:(SA{Uk)^XkU-t^Uk, 

so that BA ^XkUl® Ui and 

AdA(S) = (A, - Afe)S, (3.99) 

where the previously obtained properties of Hermitian adjoint of a simple 
tensor product of orthonormal bases were used. 

We thus have constructed a full set of eigenvectors and eigenvalues, of 
which are identically zero. This knowledge makes the practical compu- 
tation of TpiAdC) much faster and, in theory, allows of writing all terms of 
equation ( I3.93D explicitly when the characteristic polynomial is soluble. 



3.8 On some additional properties 

This section is devoted to describing how the Lyapunov exponents, and the 
system in general, behave when there are special constraints present. The 
first is simply the Hamiltonian structure, and the second is the more gen- 
eral conserved integral invariant of the flow. 

The Hamiltonian structure is usually introduced by means of a sym- 
plectic form ujs, but here, since there is already a distinguished metric 
structure, a new operator I ((J) tensor) can be used to create symplectic 
structure. We have 

ujs{X,U) = {IX,U), (3.100) 
so that for the Hamiltonian H the associated vector field is 

y = X(ttdiJ) ^ ujs{V,U) =dH{U). (3.101) 

This will also mean that there is a particular coordinate system {q,p} for 
which the coordinate basis is orthonormal and, consequently, in which the 
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connection is Euclidean, i.e. Wdq = Vdp = 0, so that 



K 

I = J2 (dPi ® - dqi J , (3. 102) 

where as before 2K = N is the dimension of the manifold. By definition the 
operator is anti-Hermitian = — J, anti-involutive 1^ = — 1 and, because 
the connection is euclidean, it is also parallel VX ~ 0. 

As follows from previous sections (ttdiJ, Z) = const and (ttdiJ, V) = 0. 
Consider next two solutions of the variational equation Zi and Z^ and the 
question of conservation 

JZ2)) = MZi, JZ2) + (Zi, J^Z2) 

= (TVz.Cttdif), JZ2) - (Zi, Vz,(ttdif)) 

= (ftdJ?, Vz.Za) - (Vz.Zi, ttdif) (3.103) 

= dH([Zi,Z2]) 

= Zi (const) — Z2 (const) = 0, 

where the appearance of the commutator follows from the zero torsion con- 
dition. The above product is thus conserved, and this means the following 
for the fundamental operator 

{FZio,TFZ2o) = (Zio,XZ2o) , 
F^IF = I, 

, 1 (3.104) 

FF^ = -FIF-^I, 

{FF^y^ = -IFIF-^ =I{FF'')I-\ 

so that the operator L defined in the previous section has the same eigenval- 
ues as its inverse, which means the Lyapunov matrix C has pairs of eigen- 
values of opposite signs {A, —A}. This also implies the flow conserves the 
phase space volume, as the sum of all such eigenvalues is zero. It should be 
stressed, that although the symplectic structure is enough to define volume 
as ?7 = LOg^, it does not give a metric structure and the respective Levi- 
Civitta connection. This is the reason for the introduction of the particular 
canonical basis {dq, dp} defined as orthonormal. 

This brings us to the next property, which is preserving the integral of 
an N-form, only this time it needs not be the volume form rj. Because the 
space of A^-forms is one dimensional the new quantity is a multiple of r/, say 
aj]. Let us assume that 

/ arj — g*s{oirj) = const, (3.105) 

Jd(s) JDa 

where D{s) is the image of some region Dq through the diffeomorphism gs. 
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Differentiating one gets 

£v{o:ri) — I d{Lv{ari)) + Lvd{ari) 



°" (3.106) 

drv(aV)ri. 

I Do 

The usual notion of divergence-free flows is just a special case when a = 1. 
To see how this definition of divergence works with the volume, and the how 
to compute with the covariant derivative, recall first an identity for the Lie 
derivative 

N 

{£vv)iEi, ...,En) = V{v{Ei, E^)) - v{Eu . . . ,[V,E,], . . . , En), 

(3.107) 

for any basis {Ei}. Together with the equality 03.71 1. and the fact that the 
torsion is zero, the above amounts to 



N 



i=l 

JL (3.108) 
dW{V)v{Ei, ...,En)=Y1 ''(-^1' ■ • • ' EU^E,V)E,, ...,E 



divF = tr(Vy), 

where the middle line is a direct consequence of the volume form being 
completely anti-symmetric and {E.,} constituting a basis. To put it shortly, 
a invariant measure exists when there exists what is called the last multi- 
plier a such that 

tr(V(aV)) = 0. (3.109) 

The special case of a = 1 (which holds for Hamiltonian systems, but not 
only), means simply that the flow conserves volume, as can be seen from 
the integral formulation above. This has a straightforward consequence on 
the Lyapunov exponents, since 

?7(Zi, ...,Zn)^ (det F)tj{Zio, Zno), (3.110) 

which holds for any A^-form. The determinant changes according to 

N 

V{t]{Zu Zn)) = J2 'K-^l^ ■ • - Zn), 

i=l 

-^(detF)?7(Zio, . . . , Zno) = tT{A)T]{Zi, .. ., Zn), (3.111) 
as 

^(deti^) =trU)detF. 
ds 

Since 2tr(£) — lndet(i^i^''') = 21n | det F\, the exponents add up to zero when 
the divergence of V vanishes. 
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Chapter 4 

Examples 



The following dynamical systems will be used to show how the Lyapunov 
exponents equation ( I3.93D can be used in practise, and also how the nor- 
mal variational equations are computed in the geometric context versus 
the algebraic one. Only the geometric considerations are in fact new, as the 
Galoisian obstructions to integrability of these systems were all analysed 
in details in the papers cited in each respective section. 

In practise, the evolution of the exponents is considered in the time t 
which is not the natural parameter, so that the results coincide with the 
standard ones. This does not change any of the formulae, as the require- 
ment that the vector field V be normalised only matters when the pro- 
jections are introduced, while the exponents are calculated for the full N 
dimensional system. The equation 03.931 1 itself is then integrated using the 
Runge-Kutta method of the fourth order. 

4.1 Arnold-Beltrami-Childress flow 

The system is given by 



We consider = and ABC =/= 0. By shifting the variables, all param- 
eters can be made positive, and as was shown in 1 11 1 the system is not inte- 
grable in the sense that there is no meromorphic first integral on the com- 
plex torus which is the system's phase space and if := C^/(2^^) < 1 
there are no real integrals either. Only one first integral is needed for in- 
tegrability, because the flow has also zero divergence and accordingly has 
a trivial last multiplier. For a three dimensional system that is enough to 
prove there must also exist two additional integrals [6|. Obviously for ^ = 
the system is separable, and solvable, and this facts allows for testing how 
the Lyapunov exponents behave. 

First, the integrable case with initial conditions a;(0) — y{0) = and 
^(0) = 0.001. Using formula 03.93D the exponents are evolved in time t, and 




(4.1) 
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can be plotted as functions of 1/i - so that the origin of the horizontal axis 
corresponds to t ^ oo. The results are presented in figures 143] and I4l2l 

A 

0.4 r 




-0.4 - 



Figure 4.1: The Lyapunov exponents for the ABC flow in the integrable case 

The maximal value of t was 10000, and the values of the exponents at 
that point are 

A = (-9.173 X 10"^ 9.099 x 10"^ 7.404 x 10"^). 

Taking now the value of /i = 1/2 to investigate a non-integrable scenario, 
and the same initial conditions, the exponents can be seen to no longer all 
be zero, there is one tending to zero and two non zero of opposite signs, as it 
to be expected for a flow with conserved volume. Their values att — 10000 
are 

A ^ (5.9 X 10"^ -5.886 x 10"^ -1.385 x lO""*), 

and their time dependence is shown in figures |4!3] and 

When it comes to the variational equation, since there are no known 
first integral (a priori), we can only go as far as projecting it with respect 
to the trajectory. As shown in [11| there is a particular solution for which 
y = 7r/4, z — x + 7r/2, so that the trajectory is described by a single equation 
in x{t) and gives as the projected VE 

r] — ^ cos(a;), 

where ^ and ry are the variations. 

In the other approach, taking the coordinate basis {dx.dydz) to be or- 
thonormal, equation ( I3.79D becomes 

^i_Z^ = {At - V<'V,)i_Z\ (4.3) 



31 



0.2 r- 




-0.2 L 



Figure 4.2: The Lyapunov exponents for the ABC flow in the integrable case 
in = 0) - a magnification of the 1/t — vicinity. 



which we want to change to involve time t, to compare it with the equations 
obtained above. Since the relation between the vectors is W ^ {W, W) V , 
the required derivations are rather cumbersome, but using the particular 
solution and the relation ds = ^J{W , W)dt, the above is finally reduced to 

= A^,i_Z'\ (4.4) 



with 



A 



sin(a;) — 2/i — sin(x) 



A = —\ 2cos(a;) -2cos(a;) | . (4.5) 
— sin(a;) 2^ sin(a:) 

This is in fact the same equation as before, when one defines £, to be ±_Z^ — 
and Tf = ±_Zy (in the cited paper the coefficient A was taken to be 1). Obvi- 
ously only two degrees of freedom are needed after the projection. 

The geometric thus gives the same "starting" equation for further inves- 
tigation. Its details - the determination of the differential Galois group can 
be found in the cited article. 



4.2 so(5) Gross-Neveu system 

The Hamiltonian reads 

H = ]:{pI+pI) - 2cos(gi) - 2cos(g2) - 4 cos(qi) 003(^2), (4.6) 
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Figure 4.3: The Lyapunov exponents for the ABC flow in a non-integrable 
case ^ = 1/2. 

which translates into the appropriate vector W 





(qi\ 


( 




Pi 


\ 


d 


92 






P2 






Pi 




-2sin(gi) 


- 4sin(9i) 


C0S(92) 






\ 


-2sin(g2) 


- 4cos((7i 


sin (92)/ 



(4.7) 



As mentioned before, this system is never integrable meromorphically 
in the Liouvifle sense as shown it Ifl2]| . As the cited paper indicates, it has 
the interesting feature of appearing regular in the original variables, and 
clearly chaotic when a complex canonical transformation q ^ iq, p ^ ~ip is 
performed. The equations presented above are those after the transforma- 
tion. In both cases the "effective" coordinates remain real, that is, when the 
initial conditions are real the variables remain real, and when they start as 
imaginary, they remain purely imaginary. 

The first set of Lyapunov exponents was obtained for the imaginary do- 
main (explicitly the above Hamiltonian) for the initial conditions of qi{0) = 
0.01, ^'2(0) = 0, pi(0) — 0.01 and P2(0) positive, determined by the condition 
H = 3. At the maximal time of 10000 the spectrum was 

A = (0.2169, -0.2169, -2.64 x 10"^ 2.64 x 10"^). 

It is clear that two of the exponents remain non zero as depicted in figure 
1431 

For the real domain, where the Hamiltonian becomes 

H = -^{pI + pj) - 2 cosli(gi ) - 2 cosh((j2 ) - 4 cosh(gi ) cosh(g2 ) , (4.8) 
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taking the initial conditions of gi(0) = 0.01, 92(0) = 0, pi{0) = 0.01 and ^2(0) 
positive such that H = —3, the exponents all tend to zero with 

A = (1.753 X 10"^ -1.753 x 10"^ -3.651 x lO"'^, 3.651 x lO"'^) 



att = 10000. The results for this case is shown in figure[ 

Let us see now how the double projection works to produce the normal 
variational equations in this case. In the algebraic approach the invariant 
plane 172 = P2 = is used to find a particular solution {qi{t) , pi{t)) , for which 
the Jacobian matrix otW is 



/ 



A = 






-6cos(gi(i)) 




(4.9) 



1 0\ 

1 



-2-4cos((?i(t)) 0/ 

Taking now only the variations in the directions 92 and p2 we get the NVE 

V = -2-4cos(gi(i)), 



(4.10) 



where ^ and ry, as before, are the variations. 

The geometric procedure is essentially the same as for the ABC flow, in 
that it consists of normalising the field W, and the appropriate equation 
03.801 1 is now 



iZ = 



(- 


-36/siii^((7i) 





-6/pl sin((7l) 


'\ 













1 




-6/pl sin((7i) 





-6/p? 





V 





-2 - 4cos((3'i) 





0/ 



Z, (4.11) 
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Figure 4.5: The Lyapunov exponents for the Gross-Neveu system in the 
imaginary domain. 



with 

r ^ 6pi(6cos(gi) ~ l)sin(gi) 
(p2-|-36sin2(gi))2 

This is, again, the same as the algebraic NVE, when two of the degrees of 
freedom corresponding to pi and qi are suppressed. Alternatively one can 
check that the two vectors with respect to which the projection takes place, 
span the pi, qi subspace, because on the trajectory 

W = (pl,0,-6sin(qi),0), 

. , , . , (4.13) 



Y = tJdff = (6sin((7i),0,pl,0). 



The above means that the next steps - checking if the NVE are soluble 
in the Liouvillian sense - is the same in both approaches. The proof that 
there are no such solutions can be found in the paper cited at the beginiing 
of this section. 



4.3 Friedmann-Robertson- Walker cosmology 

The last example is a cosmological system obtained for the FRW universe 
with a scalar field conformally coupled to gravity. It was analysed in great 
detail in |13|, and includes both integrable and non-integrable sub-cases. 
The particular Hamiltonian taken here is 

H = \{pI +pI) - \m\lql + \{kqt + Ag|). (4.14) 

When 

A = A = -m^ 
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-0.2 L 

Figure 4.6: The Lyapunov exponents for the Gross-Neveu system in the 
real domain. 

there is another first integral 

J = qiP2-q2Pi, (4.15) 

and no additional integral exists if the parameters are varied slightly. It is 
thus convenient to substitute rn^ = — Ae, so that the vector V is 





{ql\ 


( 




pi \ 


d 


q2 






p2 


dt 


pi 




-Miiqf 

-M2{ql 


+ 'zl(l + e)) 




VV 


\ 



This example is used to show, that even though the system is integrable 
only when e = 0, the Lyapunov exponents remain zero until the "perturba- 
tion" is big enough. Specifically, when e = -1, the spectrum at t ^ 10000 
is 

A = (-8.461 X 10-^ 8.461 x 10-^ -7.909 x 10-^ 7.909 x lO^^). 

Its time evolution is presented in figure 14.71 The initial conditions were 
gi(0) = 0.01, 92(0) ^ 0, pi(0) ^ 0.01 andp2(0) positive, determined by H = 
0.01. 

Changing the value of e to —1.9, with the same initial conditions gives a 
completely different picture. The exponents now read 

A = (-0.1388, 0.1388, -8.078 x 10"^ 8.078 x 10"''). 

This is to be expected for a Hamiltonian system — that the chaos becomes 
visible only for sufficient perturbation. 
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Chapter 5 

Final Matters 



The study presented here can only be considered a beginning of further 
exploration of the subject, but it can be definitely said, that the two ap- 
proaches - geometric and algebraic - can be successfully conflated to yield 
a better insight into integrability. Also, which is not to be underestimated, 
trying to use both descriptions immediately shows that some objects are 
ill-defined and some can be defined in many non-equivalent ways. 

To be more concrete, the variations or small perturbations to a given dy- 
namical system are an example of an object that reveals more when looked 
at from the geometric point of view It is fundamentally different from the 
original system itself, as it really a vector on the tangent bundle of the 
main trajectory It also turns out that it should commute with the vector 
field defining the system, as only then the transition to a nearby trajectory 
makes sense. It is also the only additional field we need to reconstruct the 
congruence of solutions if we can obtained the (first) variation as a function 
of the point on the manifold. As it is usually not the case, and we only solve 
an equation that gives the values of the variation on a particular trajectory, 
higher variations are needed. 

Here also the basic notions of differential geometry are helpful to prove 
the existence of first integrals of the higher variational equations when the 
main system has a first integral. Unfortunately it appears, that higher 
variations, when defined to agree with the algebraic definitions, are not 
coordinate independent. A fact hard to notice when analysing the equation 
only in the coordinates in which it is introduced or obtained from physical 
considerations. 

Another example of clear formulation is the Lyapunov matrix and the 
Lyapunov exponents. By definition they are constructed in a covariant way, 
but there is a price to pay for that. Namely the additional metric struc- 
ture is required. Lyapunov exponents are usually computed using the time 
parameter that is naturally present in system of physical origin, but in gen- 
eral relativistic problems or those whose formulation admits the freedom of 
time reparametrisation it is not clear which variable is the real time. And 
it is obvious that a simple exponential change of that variable could make 
positive exponents zero Il20]| . 

The calculations presented here do not require any particular choice of 
metric, so that they can be applied to any case and guarantee consistency. 
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On the other hsmd, without any particular choice it is impossible to obtain 

any results. That is why the examples included are treated as is usually 
the case - with the tacit assumption that the coordinates in which the sys- 
tem is defined are orthonormal. Until a distinguished metric structure can 
be canonically defined for dynamical systems (or at least the physical sys- 
tems), this freedom of choice will remain unresolved. 

As mentioned in the introduction there axe attempts to geometrise the 
system by finding some metric which would make the equations be the 
geodesic equations on a suitable manifold, but so far this has been done for 
a small class of systems with natural kinetic energy. It also immediately 
collides with the problem of the base space - in the case of Jacobi metric for 
example, only the configuration space is taken into account, instead of the 
whole phase space. 

This is best visible for Hamiltonian mechanics where we end up with 
analysing second order equations in the coordinates and the momenta (al- 
though also formally included in the solutions) do not play any role in the 
behaviour of neighbouring trajectories. Because the Jacobi geometrisation 
hinges heavily on the natural form of the kinetic energy, it is even impossi- 
ble to obtain an analogous picture with a space of half the dimension involv- 
ing only the momenta and suppressing the coordinates. The present work 
also shows that the symplectic structure of such systems requires some se- 
rious additional metric assumptions to speak about volume conservation. 
Even with the freedom that differential geometry gives, Hamiltonian sys- 
tems become highly structured in this context. 

Finally, among the still open problems, there is the question of studying 
more than just the Levi-Civita connection for which the results reduce to 
the algebraic ones. Introducing non-Riemannian (non-metric) connection or 
torsion, complicates the equations considerably, but has, seemingly, nothing 
to do with the question of integrability. This could hopefully give the pos- 
sibility of investigating the system on many different manifolds and in fact 
obtaining different restrictions on integrability of the same basic equations. 

Also, the algebraic tools are deeply rooted in the complex analysis of 
meromorphic functions, Riemann surfaces and analytic continuation. Thus, 
being integrable in the real sense is only understood indirectly. Here also 
lie new possibilities of extending the work to complex or Kahler manifolds, 
or developing the algebraic theory to treat the real-analytic case with more 
detail. 
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